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■ Context. The cooling histories of individual meteorites can be empirically reconstructed by using ages obtained from different ra- 

dioisotopic chronometers having distinct closure temperatures. For a given group of meteorites derived from a single parent body 
such data permit the detailed reconstruction of the cooling history of that body. Particularly suited for this purpose are H chondrites 
because (i) all of them are thought to derive from a single parent body (possibly asteroid (6) Hebe) and (ii) for several specimens 
precise radiometric ages over a wide range of closure temperatures are available. 

Aims. A thermal evolution model for the H chondrite parent body is constructed by using the cooling histories of all H chondrites for 
which at least three different precise radiometric ages are available. The thermal model thus obtained is then used to constrain some 
important basic properties of the H chondrite parent body. 

Methods. Thermal evolution models are calculated using our previously developed code (Henke et al. 2012), which incorporates 
the effects of sintering and uses new thermal conductivity data for porous materials (Rrause et al. 2011). Several key parameters 
determining the thermal evolution of the H chondrite parent body are varied together with the unknown original location of the H 
chondrites within their parent body until an optimal fit between the radiometric age data and the properties of the model is obtained. 
The fit is performed in an automated way based on an 'evolution algorithm'. Empirical data for the cooling history of H chondrites 
are taken from the literature and the thermal model is optimized for eight samples for which radiometric ages are available for at least 
I three different closure temperatures. 

Results. A set of parameters for the H chondrite parent body is found that yields excellent agreement (within error bounds) between 
the thermal evolution model and empirical data for the cooling histories of six of the examined eight H chondrites. For two of the sam- 
ples significant discrepancies exist between model and empirical data, most likely reflecting inconsistencies in the empirical cooling 
data. The new thermal model constrains the radius and formation time of the H chondrite parent body, and the initial burial depths of 
the individual H chondrites. In addition, the model provides an estimate for the average surface temperature of the body, the average 
initial porosity of the material the body accreted from, and the initial 60 Fe content of the H chondrite parent body. 

Key words. Solar system: formation, planetary systems: formation, planetary systems: protoplanetary disks 
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t^J- ' 1 . Introduction body was chosen here as a test case, because it is the only body 

for which a comprehensive set of thermochronological data is 

OO Radiometric ages for chondntic meteorites and their components available For ei ht H chondrites in particular radiometric ages 

O , provide information on the accretion timescale of chondrite par- exigt ^ were obt£dned ufji at ^ ^ difjtinct radioiso . 

CN ■ ent bodies, and on cooling paths within certain areas of these chronometers characterised by different closure tempera- 

^ . bodies. However, to utilize this age information for constraining ^ The cooli history of these samples, therefore, is well 

the internal structure, and the accretion and cooling history of constrained 
• '"J the chondrite parent bodies, the empirical cooling paths obtained 

/\ by dating chondrites must be combined with theoretical mod- The well preserved and smooth cooling histories of many 

H . els of the thermal evolution of planetesimals. Important param- H chondrites were used to argue that the H chondrite parent 

■ - - 1 eters in such thermal models include the initial abundances of body has not been subjected to catastrophic collisions during 

heat-producing short-lived radionuclides ( 26 A1 and 60 Fe), which the first * 100 Ma of its history. In that case an onion shell 

are determined by the accretion timescale, and the terminal size, structure develops in which the degree of thermal metamor- 

chemical composition and physical properties of the chondritic phism is a function of depth within the parent body; material 

planetesimals. Here we use our newly developed code for mod- of petrologic type H6 is found in the central part, while type 

elling the struc ture and thermal e volution of planetesimals of H3 material is located at the surface. Based on the assump- 

* 100 km size dHenke et al.ll2012l henceforth called paper I) to tion of an onion shell structure, a number of thermal evolution 

constrain the properties and thermal evolution of the H chondrite models were c onstructed that used increasingly more complex 



parent body, and to evaluate the most important parameters that physical input dMinster & Allegrd["l979l: iMivamoto et al.ll 1981 



determined its thermal evolution. This is done by finding the op- Bennett & McSween 1996; Akridge et al. 1998; Trieloff et al 



timal fit between a theoretical thermal model and empirical data 2003; Hevev & Sanders 2006; Sahijpal et al. 20071 iKleine et al 



for the cooling histories of H chondrites. The H chondrite parent 2008; IHarrison & Grimmll2010l) . All the previous thermal mod- 

els of the H chondrite parent body found a radius of about 

Send offprint requests to: gail@uni-heidelberg.de 100km and an accretion time of about 2 Ma after CAI forma- 
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Table 1. Key data for cooling history of selected H chondrites 



Meteorite 


Estacado 


Guarena 


Kernouve 


Richard- 


Allegan 


Nadiabondi Forest Vale Ste. Mar- 










ton 




guerite 



type 


H6 


H6 


H6 H5 H5 
Hf-W f (metal-silicate) 


H5 


H4 


H4 




radiometric age 
temperature 


4557.3 ± 1.6 
1150+75 




4557.9±1.0 4561.6±0.8 
1150+75 1100+75 






4565.6+0.6 
1070+50 


Ma 
K 


radiometric age 
temperature 






Pb-Pb g (pyroxene-olivine) 
4537.0+1.1 
1050+100 

U-Pb-Pb c (phosphates) 


4558.9+2.3 
1050+100 




4564.3+0.8 
1050+100 


Ma 
K 


radiometric age 
temperature" 


4559±2 
720±50 


4504.4+0.5 
720±50 


4522.5+2.0 4551.4+0.6 4550.2+0.7 
720+50 720+50 720+50 

Ar-Ar 6 * 1 (feldspar) 


4555.6+3.4 
720+50 


4560.9+0.7 
720+50 


4562.7+0.6 
720+50 


Ma 
K 


radiometric age 
temperature 15 


4465±5 
550±20 


4458±10 
550±20 


4499+6 4525+11 4541+11 
550+20 550+20 550+20 
Pu-fission tracks (pyroxene) 


4535+10 
550+20 


4552+8 
550+20 


4562+16 
550+20 


Ma 
K 


calculated age c 
temperature 


4465±5 
550 ± 25 


4458±10 
550 ± 25 


4499+6 4525+11 4541+11 
550 ± 25 550 ± 25 550 ± 25 
Pu-fission tracks (merrilite) 


4535+10 
550 ± 25 


4552+8 
550 ± 25 


4562+16 
550 ± 25 


Ma 
K 


calculated age c 
temperature 


4401±10 
390 ± 25 


4402±14 
390 ± 25 


4438+10 4469+14 4490+14 
390 ± 25 390 ± 25 390 ± 25 
Pu-fission tracks (pyroxene — merrillite) 


4543+20 
390 ± 25 


4544+14 
390 ± 25 


4550+17 
390 ± 25 


Ma 
K 


time interval d c 
cooling rate 


64±9 
2.5±0.4 


56±9 
2.8±0.6 


61+8 56+9 51+9 
2.6+0.5 2.9+0.5 3.2+0.6 


> -8 
(fast) 


< 8 
> 20 


12+11 
12.8+6 


Ma 
K/Ma 



Notes: 

a seelCherniak et aU ' 1 99 1 ) 
b |TrielofFetalJ<2Q03l) 

c Calculated age at 390 K from time interval between Pu-fission track (merillite, 390 K) and Pu-fission tracks at 550 K (pyroxene) compared to Ar-Ar feldspar age at 550 K 
d Time-interval for Pu- fission track cooli ng rate f rom 550-390 K, met allographic cooling rate 800-6 00 K 

e Data from U-Pb-Pb : IGopel et all 1993) . Ar-Ar: [TrTeToff et aU HW^ti . metallographic cooling rates: lTavloreraT]fT987h 

f Hf-W ages are from[KkineTniT l2008) and were re-calculated relative to the 182 Hf/ 180 Hf of the angrite D'Orbigny, which has a Pb-Pb age of / = 4563.4 ± 0.3 Ma iKleine et al|2012|) . 
g |Bouvier et 

h Recalculated for miscalibration of K decay constant (explanation see text and Trieloff et al. 2003) 



tion. The success of these models in reconstruction the empirical 
cooling histories of the H chondrites gave much support to the 
hypothesis that catastrophic collisions did not significantly dis- 
turb the initial thermal structure of the H chondrite parent body. 
It should be noted, howe ver, that some hints e xist for the oc- 
currence of impacts (see IWittmann et alj|2010i and references 
therein), but these do not appear to have resulted in catastrophic 
disruption. 

These previous thermal models, despite their success in re- 
constructing the observed empirical cooling rates for several H 
chondrites, fall short in providing insights into the physical pro- 
cesses occurring during the thermal history of ax lOOkm-sized 
planetesimals. This motivated us to use our new thermal evo- 
lution model (described in paper I) to find the optimal fit be- 
tween the empirical cooling history of the H chondrites and 
the thermal model. Our new thermal model applie s recent data 
on the thermal conductivity of granular material ( Rraus eet al.1 
I2011allbh and on the compac tion by cold isostatic pressing be- 
fore the onset of sintering (Guttl eret al.1 120091) . The modelling 
of t he sintering process is based on the same kind of the - 
ory dRao & Chakladerll 1972b as in lYoinogida & Matsuil (11984 . 
but with some improvements. This theory does not consider 
more recent approa ches for modell i ng hot pressing for techni- 
cal processes (e.g. Arzt et al .l 1 1 983b iFischmeister & Arztlll983t 
Larsson et al. 1996; Storakers et al. 1999) but appears more ap- 
propriate for the lower pressure regime relevant for asteroids. 



The treatment of heat conduction of the chondritic material and 
its sintering underpressure and at high temperature still has to be 
considered as rather crude, but at least is based on physical con- 
cepts that are successfully applied in other fields and in some 
laboratory measurements. 

We do not consider a possible melting of the body, because 
the H chondrite parent body is thought to have remained undif- 
ferentiated. Our results will show, however, that this assumption 
may not be entirely valid, because the calculated central tem- 
perature of the H chondrite parent body is only sligh tly below 
the so lidus temperature of chondritic material (1200K, iFei et al.l 
119971) . 

Our model considers more properties of the asteroidal body 
than traditional approaches, and so fitting of the model to the 
observed cooling histories is more difficult because of the much 
larger parameter space. Therefore we develop an automated 
method to perform such a fit that is based on the evolution algo- 
rithm that mimics the concepts that are thought to rule biological 
evolution in nature. This is some kind of intelligent trial-and- 
error method by which the optimum of some quality function 
can be found. Its advantage is that it does not require any kind 
of good behaviour of the quality functi on. It will be shown that 
this method (in the version described in ICharbonneaul 1 9951) can 
successfully be applied to our problem. 

For fitting the model we exclusively use meteorites for which 
at least three age determinations at different closure temperatures 
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are available. These meteorites do not only constrain the slope 
but also the curvature of an individual cooling path at a certain 
location within t he parent body. Within the framework of the 
analytic model of Miyamoto et al.] d 198 lh . such a data set for a 
single meteorite would already completely determine the radius 
and the formation time of the body. However, in view of the pre- 
sumably complex thermal evolution of meteorite parent bodies, 
a larger dataset for several meteorites is necessary for a reliable 
fit. For the H chondrites, 3 1 data from 8 meteorites (2 meteorites 
with 5 data, 3 with 4 data, 3 with 3 data) are available for this 
purpose. This allows determining all the important parameters 
of the H chondrite parent body by an optimization procedure. 

The plan of this paper is as follows: In Sect. [2] we describe 
the available empirical data. In Sect. [3]we give a brief overview 
how we calculate thermal evolution models. Section|4]describes 
the optimisation method. Section|5]gives our results and conclu- 
sions. The paper closes with some final remarks in Sect. [6] 

2. Data for cooling histories of H chondrites 

Most ordinary chondrites were affected by thermal metamor- 
phism in their parent bodies and exhibit a range of meta- 
morphic condi tions from type 3 (unequilibrated) to type 6 
(equilibrated) dDoddl 119691) . This range in metamorphic con- 
ditions probably reflects different heating and cooling histo- 
ries, presumab ly caused by varying burial depths within the 
parent bod y (Miya moto et al.lll98UfMcSween & Patchen|[l989t 
ITrieloff et alj|2Q03h iKleine et al.ll2008l) . The heat source for the 
thermal metamorphism of ordinary chondrites most probably 
was heating by decay of 26 A1 and to a lesser extent also 60 Fe. 

Ordinary chondrites and their components h ave been dated 
by a variety of methods dNvquist et al.l [2009). Al-Mg ages 
for chondrules from L and LL chondrites indicate that chon- 
drule formation occurred approximately 2Myr after forma- 
tion of Ca-Al-rich inclusions (CAIs commonl y regarded as the 
oldest objects formed in the solar system) (iKita et"al] [2000: 
iRudra swami & Goswarnl 120071: Ikudraswa mi et alf 20081) . The 
formation age of chondrules from H chondrites is less well 
known because no primitive H chondrite has yet been investi- 
gated with the Al-Mg system. However, a Hf-W age for the H4 
chondrite Ste. Marguerite appears to date chondrule formatio n 
at 1.7 ± 0.7 Myr after CAI formation (IKleine et al.ll2008l 12009b . 
about contemporaneously to chondrule formation in L and LL 
chondrites. 

The ages for chondrules from ordinary chondri tes mainly 
constrain the timescale of parent body accretion (lAlexanderl 
12005b Kleine et al . 2008). In contrast, ages for equilibrated chon- 
drites can be used to reconstruct the cooling history of ordinary 
chondrite parent bodies, by applying chronometers having dif- 
ferent closure temperatures to a set of equilibrated chondrites of 
different metamorphic grade. By far the most comprehensive set 
of such age data exists for H chondrites (see Table [TJ, providing 
temporal information for cooling over almost the entire range 
of temperatures prevailing from the metamorphic peak down 
to ambient temperature. The closure temperatures for diffusive 
exchange of parent an d daughter elements among the differ- 
ent minerals in a rock dDodsonlll973t iGangulv & Tirond l2001) 
for radiometric methods applied to H chondrites range from 
» 1150K (Hf-W) to * 550 K (Ar-Ar). Additional information on 
the cooling history is provided by 244 Pu fission tracks in orthopy- 
roxene (corresponding to cooling below w 55 K) and merrillite 
(corresponding to cooling below w 390 K) dPellas et al.lll997t 
ITrieloff et al.l l2003). and from metallographic cooling rates (cor- 
responding to cooling between * 800 and m 600 K), which em- 



Table 2. Fixed parameters of the H chondrite parent body mod- 
els (taken from paper I). 



quantity 


value unit 


bulk density pt, 


3.78 gem- 3 


mass fraction X A] 


9.10 x irr 3 


mass fraction Zp e 


2.93 x 10-' 


26 A1/ 27 A1 ratio (at CAI formation) 


5.1 x 1(T 5 



ploy Fe-Ni diffusion pr ofiles in metal grains consisting of ka- 
macit e and taenite (e.g., Her pfer et al.lll994l:lHopfe & Goldstein! 
120011) . 

The age and thermochronological data available for H chon- 
drites are summarised in Table [1] along with information on the 
closure temperatures of the different methods. These age data 
permit a detailed reconstruction of the cooling history of H chon- 
drites, and form the basis for numerically simulating the accre- 
tion and cooling history of the H chondrite parent body (Sect.|5]l. 

For calculating heating rates the difference between the for- 
mation time of CAIs and the ages from TableQ]has to be formed. 
This sets the zero-point of our time-scale and the instant for 
which the initial abundance of 26 A1 is prescribed . As reference 
time w e take the value of 4 568.5+0.5 Ma from iBouvier et ail 
(2007). 



3. Thermal evolution models 

The temperature history of an undifferentiated meteoritic parent 
body is assumed to be determined by essentially three processes: 
(i) Transient heating by decay of short-lived radioactive nuclei 
26 A1 and 60 Fe and long lasting heating by long-lived radioac- 
tives ( 40 K, 232 Th, 235 U, 238 U), (ii) transport of heat to the surface 
by heat conduction, and (iii) energy exchange with the environ- 
ment by emission and absorption of radiation energy. The heat 
transport in the body is governed by the fact that the material is 
porous, that heat condu ctivity of chondritic material strongly de- 
pends on porosity (e.g. lKraus e et al. 201 lb), and that the poros- 
ity changes by sintering as the body is heated up. 

For the formation of the bodies we apply the "instanta- 
neous formation" approximation, i.e., it is assumed that the bod- 
ies acquired most of their mass during a period short com- 
pared to the half-life T\j% - 7.2 x 10 s a of the main heat- 
ing source, 26 A1. Present planetary formation scenarios suggest 
such rapid mass acquisition of 100 km sized and bigger b odies 
(e.g.lWeidensc hilling & Cuzzil2006tlNagasawa et al.l2007|) . The 
rapid growth is approximated by the assumption that the bodies 
came into life at some instant ff orm and have constant mass over 
their subsequent evolution. 

For reasons of computational economics the shape and inter- 
nal structure of the body is assumed to be spherically symmetric, 
though this is certainly a strong simplification. In detail, the fol- 
lowing equations are solved: 

1 . The equation of hydrostatic equilibrium 

dp GM r 

1~ = 2"^' f 1 ) 

d r r z 



for the pressure p, where 

M r = 4n f dr' r' 2 g(r') . 
Jo 

and 

Q = QbD . 



(2) 



(3) 
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The density q\, is the mass-density of the consolidated material 
and D the fraction of volume filled with matter. The void fraction 
or porosity of the material is <f> = 1 — D, 

2. The heat conduction equation for the evolution of the tem- 
perature 



dr 



1 d 



P A 1 



c v 1 -— — r q r — +h h v r F r 



d r gr 2 d r 



p At Q 



(4) 



where c v ist the specific heat of material per unit mass, h the 
specific heating rate by decay of short living radioactives, v r F r 
the local gravitational acceleration and v r the shrinking velocity 
during sintering. The heat current is 



8T 
or 



(5) 



with the heat-conduction coefficient K. 

The heat conductivity used in our model depends on the 
porosity <p = 1 - D of the material. We use the analytic fit 



K(D) 



Kh ^ e -4(l-D)/0i + e 4(a-(l-D)/fc) ^ 



1/4 



(6) 



to experimental data discussed in iKrause etal ] (1201 lah and in 
paper I, where <pi = 0.08, <p2 = 0. 167, and a = —1.1 are constants 
given in paper iFlThe pre-factor corresponds to the heat con- 
ductivity of the bulk material (i.e., at D — 1). The average value 
of K\, is determined by e xtrapolating data obtained from a couple 
of ordinary chondrites ( Yomogida & Matsui 1983]) by means of 
Eq. © to zero porosity. Its value is K — 4 W m s _1 , but there is 
significant scatter. 

3. An equation for the sintering of the initially porous mate- 
rial under the influence of pressure 



3D 

— = F(D,p,T), 
at 



(7) 



where the rhs. is determined b y solving a set of equation s for 
the specific sintering model of iKakar & Chakladerl (1 19671) and 



iRao & Chakladeri (119721) . 

The internal structure and thermal evolution of parent bodies 
of planetesimals is calculated by the numerical method described 
in paper I. 

The radius of the body is not fixed but changes during the 
evolution by sintering of the initially porous material. If the body 
would be completely sintered the final radius would be R. This 
radius is related to the mass by M = (4^/3)pb^ 3 with p\, be- 
ing the bulk density of the completely consolidated material. 
Initially the material is porous and the radius is bigger because 
the density of the porous material p = Dp\, is lower. The initial 
porosity is <f> SI f. This equals also the porosity of the outermost 
surface material which is too cool and under too low pressure 
for sintering. 

The problem depends on a number of parameters, the most 
important of which are the bulk density pb of the chondritic ma- 
terial, the mass M of the body or its radius R, the formation time 
fform (measured relative to the formation time of CAIs), the heat 
conduction properties of the material (.Kb), the mass fractions 
Xai and Xp e of Al and Fe, respectively, in the chondritic mate- 
rial and the quantities of short-lived radioactives 26 A1/ 27 A1 and 
60 Fe/ 56 Fe the body is endowed with, the porosity SI f of the sur- 
face material, and the surface temperature T SI f. 

The surface temperature T SI f depends during the initial phase 
of evolution on the properties of the ambient dusty accretion 



1 In paper I the constant a was unfortunately given without the minus 
sign 




Fig. 1. Schematic representation of the temperature evolution at 
two different depths (full lines) and at the centre (dotted line) of 
a planetesimal of the 100 km size class. The horizontal dashed 
lines correspond to fixed temperatures equal to the closure tem- 
peratures in Table [TJ The circles indicate the instants at which at 
a given depth the temperature drops below a given temperature. 
These instants correspond to the closure ages of the different 
thermochro no meters used for meteorites. 



disk, after disk dispersal on the distance to the proto-sun. Since 
model properties depend only weakly on r sr f, this quantity is 
held fixed over the whole evolution period in this study. 

Some of the parameters are fixed by the properties of the 
H chondrite material and general properties of the solar system. 
They are given in Tableland their values are chosen as in pa- 
per I. The remaining parameters M or R, ?f olm , T sr f, <p ST f, and K\, 
are special properties of individual bodies and are varied in our 
model calculations. 

Since there seem to be doubts that the 6l) Fe/ 56 Fe ratio 
(at time of CAI formation) is the same for all bodies (see 
Ouitt e et al.ll2010h we prefer to consider also this quantity as a 
free parameter (details are given in Sect. 15. U . 



4. The fitting procedure 

4.1. Comparing thermochronological data with asteroid 
models 

The data on the cooling history of meteorites described in Sect. [2] 
associate with a meteorite a number of temperature-time data 
pairs, corresponding to the instants where the temperature at the 
pristine location of the meteorite inside of its parent body fell 
below that temperature. The set of data pairs obtained for a me- 
teorite corresponds to a set of intersection points between the 
temperature evolution curve of some mass-element with the lines 
of constant temperature on the descending part of the evolution 
curve. This requires, of course, that the temperature evolution at 
the location of the later meteorite proceeds without substantial 
disturbance by external influences (impacts, for instance). 

Figure[TJdemonstrates this for the four closure temperatures 
( 1 150 K, 720 K, 550 K, 390 K) of Tabled] The cooling curves are 
shown as full lines, the lines of constant temperature as dashed 
horizontal lines. The circles on a cooling curve correspond to 
the data for a hypothetical meteorite. Since the average slope of 
the curve connecting two neighbouring data points significantly 
varies for curves corresponding to different depth's, it is obvi- 
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Table 3. Parameter space of variable parameters for the mete- 
oritic parent bodies 



quantity 


value 


unit 


Formation time tf onn 


1.5-3 


Ma 


Radius R 


0-200 


km 


60 Fe/ 56 Fe ratio 


0-2.0x 10- fi 




Initial/surface temperature T SI { 


150-400 


K 


Surface porosity (f> slf 


0.2-0.8 




bulk heat conductivity A'b 


1.0-5.0 


Wm-'r 1 



ously possible to reconstruct for a given model from only two 
data points the appropriate depth for which the temperature evo- 
lution curve runs through just these two data points. 

If the properties of the body are not known and only two data 
points are available for a meteorite, such a fit would be possible 
for a whole range of model parameters. This does not allow any 
conclusion on the properties of the body. If, however, more than 
two data points are available for a meteorite this restricts the pos- 
sible models since the overall shape of the temperature evolution 
curves depends on the model properties. The average slope of 
that part of the curve that connects neighbouring data points on 
the curve varies for different such pairs with model properties 
in a different way. This allows to reconstruct the properties of 
the body from cooling data on meteorites by looking for a set 
of model parameters for which for each of the meteorites an ap- 
propriate depth can be found such that the cooling curve runs 
through all data points. 

Since for real meteorites, experimentally determined closure 
temperatures and ages are subject to considerable errors (cf. 
Table[TJi, data sets for an as big as possible number of meteorites 
(with at least three data points) are required and the task is to find 
a set of parameters that describe the properties of the parent body 
such that for each of the meteorites a depth below surface can be 
found such that the (somehow defined) minimum distance of all 
of the data points of a meteorite from the corresponding tem- 
perature evolution curve becomes as small as possible; and this 
should hold for all available meteorites. Hence one has to solve 
a minimisation problem, where some quality function has to be 
minimised with respect to a set of data and a set of parameters. 

bf In the past, this task has been solved mainly by trial and 
error. Formation times ff orm and parent body radii R for some 
bodies have been de t ermin e d this way ([Mi yamoto et al. 198 1 ; 
Bennett & McSween 1996; Akridge et al. 1998; Trieloffetal 



2003 



2008 



Hevev & S ande rsll^OOalS ahii pal et ai.l l2Q07b iKleine et al 



Harrison & G rimm 2010). However, the essential proper- 



ties of an asteroidal body that determine its thermal evolution, 
are not merely R and ff orm but also some additional parameters, 
see Table [3] To these, one has to add the unknown depth's for 
each meteorite as additional parameters of the optimisation prob- 
lem. Here we aim to solve the problem by a general numerical 
algorithm that allows to treat the more extended parameter space 
(and possibly additional parameters in future calculations) and 
large data sets. 

4.2. The data set 

For reconstructing the parent body of a group of meteorites that 
are thought to derive from the same asteroidal body, we use 
all specimens for which at least three thermochronological data 
points are available. For the parent body of, e.g., the H chon- 
drites, the eight meteorites given in Table [TJhave sufficiently ac- 



curate determined data and the required minimum number of 
data points. Three of them have even four data points, the other 
five have three of them. 

The burial depth of each of the J meteorites, enumerated by 
an index j (= 1, . . . , J) , in the parent body is characterised by 
the J mass-coordinates Mj in the theoretical model of the parent 
body, because the model calculation uses the mass-coordinate as 
independent variable (see paper I). Then we have a list of / data 
quintuples 



TUj = (Tj, <T T ,i, tj, <T,J, jj) , (/=!,...,/) 



(8) 



with Tj being the closure temperature, f, the corresponding age, 
(jj j and o~ i j the errors of Tj and f,, respectively, and j, a pointer 
to the appropriate member in the list (Mj, j e { 1 , . . . , J}) of mass- 
coordinates corresponding to the different meteorites. 

The mass-coordinates Mj are not known a priori. They have 
to be determined simultaneously with the parameter set given in 
Table[3]by the optimisation process. 

4.3. The quality function 

The quality function to be optimised can be defined in differ- 
ent ways. We construct in this paper one in the following way: 
For each of the mass-coordinates Mj it is known which of the 
data sets cr, are associated with this Mj. During a model run, 
the instants when the temperature drops below the closing tem- 
peratures belonging to the data sets m-, associated with j, are 
recorded. This defines for each i an age f mo d,i on the temperature 
evolution curve j where it passes through the prescribed tem- 
perature Tj. Analogously, the temperature r mo d,, is recorded at 
instant f,- for each of the data points. 

As quality function for the optimisation process we choose 
the least squares deviation of the model age f mo d,, from the em- 
pirically determined age f,, and of the model temperature T mo ij 
from the closure temperature T, : 



;=1 V 



t{ Anod,< 

°~t,i 



I — I = * 



0~T,i 



(9) 



This x~ is to be minimised. The parameters (R, tf orm , 60 Fe/ 56 Fe, 
T sr f, <p sr f, Kb, Mj [j — 1, . . . , J)) are then determined such that 
each cooling curve corresponding to the burial depth of a certain 
meteorite in the model passes as close as possible to the em- 
pirically determined age-closure temperature pairs for all data 
points of that meteorite. 

4.4. Choice of the optimisation algorithm 

The basic problem is to determine the absolute minimum of the 
quality function O(Tir) and possibly all other deep local potential 
minima with comparable depth. As is well known, such an op- 
timisation problem with a large number of variables cannot be 
solved by analytical methods and the task of solving the prob- 
lem by numerical methods becomes extremely difficult and time 
consuming if the surfaces <S>(m) = const have a complicated ge- 
ometric structure, which is expected to be the case for our prob- 
lem. No general method is known which allows to definitely 

1 . find all local extremal values of an optimisation problem and 
to 

2. decide if the real optimum indeed is contained within the set 
of extremal values found by the applied method. 
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Methods have been developed which allow one to find an in a 
certain sense best result, which means that for test cases with a 
known absolute minimum the algorithm always finds this mini- 
mum or at least a local minimum that is not much worse than the 
absolute one. 

One method, that has successfully been applied to several 
difficult optimisation problems like that of the tra velling sales- 
man is the "evolution algorithm" developed by iRechenbergl 
(Il973l[l989l) . A detailed description of how this method can be 
applied to p roblems in astronom y and astrophysics and a code 
are given in ICharbonneaul (1 1 9951) . We use this method to recon- 
struct the parent body's properties of a suite of meteorites. 

The basic idea is to imitate the concepts that are thought to 
rule the evolution of species in an algorithm to find the extremal 
values of some quality function. Evolution theory describes the 
change of the shape of creatures in a population during many 
generations due to the variability of the population and the de- 
pendence of the amount of the creature's offspring on their con- 
formity to the environment. The biologic term genotype means 
the whole heritable attributes of a creature, while phenotype de- 
notes its appearance. During the development of the creature the 
phenotype forms according to the instructions of the genotype. 
The position in the genotype, where a certain property (e.g. eye 
colour) of the individual is encoded, is called gene. The value of 
the property (e.g. blue) is called allele. The amount of possible 
offspring (the reproduction outcome) depends on, how advanta- 
geous its alleles are. The reproduction outcome determines, how 
numerously certain alleles in a population are inherited com- 
pared to other respective alleles. Which processes are advanta- 
geous is determined by the environment. Thus after some gen- 
erations advantageous alleles will accumulate in the population, 
while less advantageous ones will be depleted, such that after 
a certain number of generations only those alleles, the creature 
benefits most from, remain in the population in a noteworthy 
amount. This leads to a degeneration of the gene pool. But ad- 
ditionally mutations can cause random change of the genotype 
in a certain allele, which can cause change of the affected gene 
in a significant way and can create a completely different pheno- 
type. Thus mutation adds novel alleles to the gene pool. These 
can be superior to the established alleles (even if they are de- 
structive in most cases). So the mutations keep the variability 
of a population, and allow the population to adapt to changing 
environments. 

An algorithm for finding parameters that works like evolu- 
tion due to inheritance and natural selection is superior in finding 
of optimal parameters in ugly problems with many local min- 
ima compared to methods that, e.g., calculate in some way the 
derivative of the quality function and then go into the direction of 
downward slope and repeat calculating the derivative and so on. 
Those methods readily find the next local minimum which, how- 
ever, might not be the absolute one. In the evolution algorithm 
mutations can produce extensive changes of the parameters to 
the parents' ones. This can cause big jumps in the parameter 
space. Thus the population will not remain in a local minimum, 
but test remote regions again and again and perhaps find a deeper 
local minimum that is possibly a total one. 

4.5. General implementation of an evolution algorithm 

An evolution algorithm can be applied to every problem in which 
the variables to be optimised can be encoded to form a string 
of numbers (chromosome). This string then is the genotype. It 
is assumed in the following that for every individual the model 
is calculated and compared to the data by calculating a quality 



function, like Eq. (|9), e.g., and ranked by this to other individ- 
uals. After determining which individuals are to get off-springs, 
their genotypes are cut at the same random digit(s) and the re- 
lated parts are exchanged, so that, e.g., the initial part of parent 1 
is connected to the end part of parent 2. The so generated string 
then represents the genotype of the new offspring individual. 
Mutation can be performed by giving every digit of the string 
of the genotype a certain possibility to take a random value. 

An implementation of an evolution algorithm can be as fol- 
lows: 

1 . A random initial population is constructed with normally no 
reference to the problem. 

2. The fitness function of each individual is calculated by a user 
given algorithm which determines how well the individuals 
match the problem. 

3. A new generation is constructed by breeding of selected in- 
dividuals of the old population. There the algorithm creates 
new parameters by exchange of information between the pa- 
rameters of the two parent individuals. Also the mutation 
takes place in this step. In most cases the number of indi- 
viduals in each generation will remain the same. 

4. The old population is replaced by the new one. 

5. Then back to 2. 

4.6. Application to our problem 

The particular impl ementation of a gene tic algorithm that we use 
is that described in ICharbonneaul (I 1 995b . It uses a random num- 
ber generator that generates numbers between and 1 to derive a 
set of random parameter values for formation time, radius and all 
other parameters of the model. Each time the genetic algorithm 
requires a new value of O for some individual it provides a set of 
K random numbers < f * < 1 (k — 1 , . . . , K) from which values 
for the corresponding physical parametres x^ of the model are 
calculated as 

%k — Xk.mm (^ft.max -^min) & ■ (10) 

Here jc* im i n and Jt^max are the minimum and maximum allowed 
values for xi according to Table [3] For this set of parameter val- 
ues Xk an evolution model is calculated and from its results the 
corresponding value of the quality function <£. 

The burial depth's of the meteorites are treated as part of the 
parameter set of the model. Random values are assigned to the 
burial depths of each of the distinct H-chondrites. In this way 
for each of the individuals of the evolution algorithm, a certain 
depth is ascribed to every chondrite. Technically this is done by 
assigning to each meteorite a layer characterized by an index i 
and a position factor y that gives its relative position inside of the 
layer with respect to the edges of the layer, r, and r,_i , so that the 
real position of the meteorite is 

r me t = n-i + y ■ (n - . (11) 

Because the sintering behaviour is treated as the same every- 
where within a single layer the factor y will not change during 
the thermal evolution. The temperature at the position of the me- 
teorite is calculated by 

T met = T i - l +y-(T i -T i - 1 ) (12) 

and it is tested how good the temperature evolution of that r met 
matches the empirical data of the corresponding meteorite by 
means of Eq. (0. 
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Table 4. Some properties of asteroid (6) Hebe. 



(J LldllLlLy 




H i m pn c i on 
U.1111C11MU11 




average radius R 


93 


km 


a 


total mass M 


1.28 1 22 


g 


b 


average mass-density g 


3.81 


g cirT 3 


b 


radius of orbit a 


2.462 


AU 


c 


excentricity e 


0.23 




c 


inclination i 


14.8° 




c 



Sources: fa OTorppa et alj t2003) . (b^ lBaer et al J t20TH) . (c)|Gafrev & Gilbert GUI) 



According to Charb onneaul dl9 95 ) a fit achieved by a genetic 
algorithm is moderately good if I—K - x 1 , where / is the number 
of data points to be fitted and K the number of the optimised 
parameters. For example, for the six parameters given in Table[3j 
plus the eight layers of the eight meteorites (K — 6 + 8 = 14), 
and the 28 data points, one gets a desired value of^ 2 < 14 for 
an acceptable fit. 



5. Results for the parent body of H-chondrites 

5.1. Parameters 

Besides the physical parameters of the problem, that are not sub- 
ject of the fitting procedure, we have a number of parameters that 
describe properties of the body and external quantities. These are 
given in Table [3] including the range of variation admitted dur- 
ing the optimisation process. The following parameters require 
some comments: 

R adius: There are a num ber of reasons (iGaffev & Gilber j 

119981: iMorbidelli et all 1 19941) why asteroid (6) Hebe is thought 
to be the parent body of the H chondrites and the HE iron me- 
teorites. Some general characteristics of (6) Hebe are listed in 
Table [4] If (6) Hebe is only a fragment of the original H chon- 
drite parent body, then the size of (6)Hebe provides a firm lower 
limit for the size of H chon drite parent body. Hebe is a some- 
what irregular shaped body ( Torp pa et al.l2003l) with dimensions 
205 x 185 x 170 km corresponding to an average diameter of 
186 km. Its high mass-density means that it is not a rubble pile 
as most other asteroids. It did not suffer catastrophic collisions 
that disrupted the whole body and where the debris later re- 
assembled to the present body. Nevertheless, a surface layer of 
unknown thickness has certainly been eroded over time by colli- 
sions with other bodies and the present size of Hebe defines only 
a lower limit to the size 4.56 Ga ago. But this lower limit should 
be observed in the model construction. 

60 Fe/ 56 Fe-raf/o: The decay of 60 Fe may have contributed 
to the heating of the body, in particular sin ce the recently re- 
determined half-life of r 1/2 = 2.62 ± 0.04Ma dRugel et al.f2 009) 
significantly exceeds the half-life of 26 A1 of T1/2 = 0.72 Ma. 
Though there is less 60 Fe than 26 A1, the heating period can be 
significantly prolonged in models with 60 Fe compared to models 
with only 26 A1 included. 

It has been questioned (cf. Dauphas et al. 2008^ lOuitte et"aT1 

l2010l) that 60 Fe was homogeneously distributed throughout the 
solar system. If true, this would mean that the initial 60 Fe abun- 
dance in the parent bod y is not known a priori. On the other 
hand, a recent study by iTang & Dauphasl (1201 2l) argued for a 
homogeneous distribution of 60 Fe on a relatively low level of 
1 .08+0.21 x 10 -8 for the initial 60 Fe/ 56 Fe abundance ratio at time 
of CAI formation. This is much lower than values previously re- 



ported, although the authors do no comment on the credibility of 
earlier studies. 

At such a low abundance heating, by 60 Fe is nearly negli- 
gible. Because of the presently unclear situation with respect to 
which value should be taken for the 60 Fe/ 56 Fe abundance ra- 
tio we treat this as an additional free parameter for the optimi- 
sation process. The range of 60 Fe/ 56 Fe ratios considered (see 
Table [3]) is betwe en no 60 Fe and the maximum reported value 
dBirck&Lugmairj|1988l) . 

Surface temperature: This is determined by the unknown ra- 
dius of the orbit, its eccentricity, by the rotation of the body, the 
surface albedo, and, during the first few million years, by the 
properties of the ambient accretion disc. All these factors are not 
known precisely. Therefore we think it best to include the surface 
temperature T st f as a free parameter in the optimisation process. 
The resulting temperature r srt then represents some kind of aver- 
age temperature of the surface layers during the evolution period 
where the body first heats-up and then finally cools down below 
the lowest of the closure temperatures. 

The resulting value for r srt has no bearing to the present sur- 
face temperature of a possible successor body, because in the 
meantime the body may have changed its properties by catas- 
trophic collisions and it may have changed to some extent its 
orbit during the evolution of the asteroid belt. Surface temper- 
atures in the present-day asteroid belt therefore cannot be con- 
sidered as representative for the surface temperature during the 
phase where the essential metamorphic processes acted in the 
body, as it is assumed in some thermal evolution models. 

The range of surface temperatures (see Table[3} is limited by 
the fact that T SY f should be higher than the temperature where the 
body would contain ice (« 150 K) and lower than the lowest of 
the closing temperatures (« 400 K). 

Surface porosity: During the (short) main growth phase of 
the body the structure of the surface layer is determined by the 
properties of the smaller planetesimals that form the body under 
consideration. Later the surface layer is gardened by continued 
impacts that produce a regolith surface layer for which theoret- 
ical calculations indicate a poro sity in the range 20% to 30% 
(Rich ardsonll2009l : IWarrenll201 ll) . Presently, too little is known 
on these processes as that they can be included in our model cal- 
culation. Therefore we include the porosity of the surface layer, 
sr f, as a free parameter in the optimisation procedure. The re- 
sulting value then is some sort of average value of sr f over the 
period where the metamorphic processes where active. 

5.2. Data 

The general optimisation method of the evolution algorithm is 
applied to the data for the set of H chondrites given in Table Q] 
For reasons that will become clear from the model results shown 
in the subsequent figures and are discussed later the two highest 
closing temperature data of Ste. Marguerite and the plutonium 
fission track datum of Nadiabondi are not used. Then we have 
totally 28 data points that have to be fitted. The problem depends 
on the 5 parameters given in Table. [3] and the eight unknown 
burial depths of the meteorites. 

The range of surface porosities (see Table [3]) should in prin- 
ciple be restricted by the requirement, that the porosity should 
be higher than that of the m ost porous chondrites (about 20%, 
cf. IConsolmag no et alj |2008). During the model calculations it 
was found that for an optimum fit values not much bigger than 
this minimum value are required and, therefore, no a priori re- 
striction was imposed on <f> SI f during the final optimisation runs. 
The optimal set of parameters, indeed, satisfies the assumption. 
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Fig. 2. Results for the individual meteorites used in the model fit. The values of the free parameters and the burial depth's of the 
meteorites found for the optimised model are given in Table [5] The full line shows the temperature evolution at the burial depth's 
of the meteorites in the asteroid, as determined by the evolution algorithm. The dashed line shows the temperature evolution at the 
centre and the dotted line is the surface temperature. The square boxes indicate the individual data points and the error bars the 
uncertainty of closing temperature and age determinations. 



5.3. Model with fixed heat conductivity 

An optimisation was run for the 5+8=13 parameter and the 28 
given data points. The model used 250 mass-shells. The optimi- 
sation was run for about 5 500 generations with hundred indi- 
viduals per generation, i.e., totally a number of about 550000 
complete thermal evolution models were calculated for a com- 
plete optimisation run. This used about seven days on a normal 
desk-top computer equipped with a dual-core pentium processor 
operating with 2.6 GHz. 

The finally accepted optimum fit had a value ofy 2 = 5.881 in 
generation 2987. This value for^ 2 can be considered as excellent 
in view of the fact that a good fit would require only^- 2 < 15 (cf. 
Sect. I4.61 l. The model parameters of the final optimum fit for the 
properties of the body and the burial depths of the meteoritic 
material are given in Table [5] 

The central temperature at maximum, T c — 1 187 K, that re- 
sults for the set o f optimum param eters is below the Fe-FeS eu- 
tectic at 1220K dFei et all [l 9971) . For the optimum model the 
temperature in the centre therefore is below the border tempera- 
ture where differentiation would just become possible. In so far 
the model is compatible with the basic assumption of the model 
calculation that melting can be neglected. 



The cooling curves of the meteoritic materials at the result- 
ing burial depths and the empirical age-temperature data points 
given by the radiometric ages are shown in Fig. [2] for the op- 
timised model. The solution of the optimisation problem found 
here meets the data for the cooling history of the meteorites of 
Table[T]within their lcr-error-boxes in all eight cases, though in a 
few cases only marginally. The fit between model and empirical 
data can be considered as excellent for the three meteorites with 
four data points: Kernouve, Richardton, and Estacado. Also for 
three of the meteorites with three data points, Guarena, Allegan, 
and Forest Vale, the model fits the data-points rather well. Only 
for Nadiabondi and Ste. Marguerite the fit is of low quality. 
These two cases are discussed later. 

The maximum temperature, r max , reached during the ther- 
mal history at each of the different depths are given in 
Table [5] For comparison, the table also gives the meta- 
morphic temperature of the differ ent petrologic types, r met 
dSlater-Revnolds & McSweenl 120051) . The maximum tempera- 
tures of the different meteorites at their burial depth, T mllx , fit 
reasonably well with the requirement r max > r met . 

The radius of 140 km required for the optimum model is sig- 
nificantly larger than the present day radius of (6) Hebe. If this 
asteroid is really the parent body of the H-chondrites this would 
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Fig. 3. Results for a model fit where the heat conductivity at zero porosity, K\,, is included in the optimisation process. The values 
of the free parameters and the burial depth's of the meteorites found for the optimised model are given in Table [6] For meaning of 
symbols see caption to Fig. [2] 



Table 5. Parameters and burial depths of the meteorites for the 
optimised H chondrite parent body model 



Quantity 




value 


unit 




Radius R 




140.0 


km 




Formation time t folm 


1.927 


Ma 




Heat conductivity K b 


4.00 


Wm-'K" 1 




Surface temperature T SI f 


241.2 


K 




60 Fe/ 56 Fe ratio 




0.0 






Surface porosity S 


r 


20.0% 






Maximum temperature T c 


1 187 


K 




Porous outer layer 




1.166 


km 




Meteorite 


type 


depth (km) 


Tmax (K) 


Tract (K) 


Estacado 


H6 


61.8 


1174 


S3 1170 


Guarena 


H6 


60.5 


1174 


ss 1170 


Kernouve 


H6 


48.5 


1169 


ss 1170 


Richardton 


H5 


18.8 


1090 


S3 1020 


Allegan 


H5 


17.6 


1079 


« 1020 


Nadiabondi 


H5 


15.9 


1065 


S3 1020 


Forest Vale 


H4 


5.60 


886 


S3 970 


Ste. Marguerite 


H4 


1.60 


753 


«970 



mean, that a layer with an average thickness of nearly 50 km has 



been eroded over time from the surface of (6) Hebe. This is more 
than the depth of the layer where in this model the petrologic 
type changes from H5 to H6. The change from initial to present 
radius amounts to 33% or corresponds to a volume change of 
71%, i.e., the body would have lost most of its initial mass and 
all H3 and H4 chondrite material should have been removed, 
which is inconsistent with meteoritic record. Probably our op- 
timised radius is too big and the most likely reason for this is 
that the treatment of heat conduction in the porous medium and 
of the sintering process is not yet realistic enough. The other 
possibility is that H-chondrites do not come from (6) Hebe, al- 
though ^fhe_argjim^nts_fevouring (6) Hebe as parent body are 
strong dGaffev &Gilbertlll998h . 

The thickness of the remaining porous outer layer is also 
given in Table [5] As lower boundary of this layer we define the 
radius where (p(r) = 10%. At this degree of porosity usually 
the connections between pores start to close and the pore struc- 
ture changes from a network of interconnected pores to isolated 
voids. 

The low value of the surface porosity sr f = 0.20 may either 
result from the fact, that planetesimals of the 100 km size-class 
gain most of their mass by collisions with smaller planetesimals 
that have already diameters of about 10 km and more and, there- 
fore, are already largely compacted (see Fig. 10 of paper I) if 
they are incorporated into the growing larger body. A second 



9 



Stephan Henke et al.: Thermal history modeling of the H chondrite parent body 




Fig. 4. Results for a model fit where the radius is set equal to the observed radius of (6) Hebe. The burial depth's of the meteorites 
found for the optimised model are given in Table[7] For meaning of symbols see caption to Fig. [2] 



Table 6. Parameter for the model with optimised Kq, burial 
depths of the meteorites, and maximum temperature during their 
thermal history. 



Table 7. Parameter for the model assuming a radius equal to 
that of (6) Hebe, burial depths of the meteorites, and maximum 
temperature during their thermal history. 



Quantity 




value 


unit 


Quantity 




value 


unit 


Radius 


R 


123.6 


km 


Radius 


R 


93.0 


km 


Formation time 


fform 


1.880 


Ma 


Formation time 


^form 


2.239 


Ma 


Heat conductivity 


K h 


2.65 




Heat conductivity 




2.92 


Wnr'KT 1 


Surface temperature 




225.0 


K 


Surface temperature 


r srf 


306.5 


K 


60 Fe/ 56 Fe ratio 




2.9 x 10~ 7 




60 Fe/ 56 Fe ratio 




2.0 x 10~ 7 




Surface porosity 


0srf 


29.4% 




Surface porosity 


0srf 


34.0% 




Maximum temperature 


T a 


1234 


K 


Maximum temperature 


T 

1 c 


1 154 


K 


Porous outer layer 




0.396 


km 


Porous outer layer 




0.525 


km 


Meteorite 


type 


depth [km] 


Tmax [K] 


Meteorite 


type 


depth [km] 


T mia [K] 


Estacado 


H6 


43.1 


1212 


Estacado 


H6 


69.7 


1 149 


Guarena 


H6 


43.2 


1212 


Guarena 


H6 


70.9 


1 150 


Kernouve 


H6 


34.0 


1203 


Kernouve 


H6 


47.9 


1 133 


Richardton 


H5 


13.7 


1 111 


Richardton 


H5 


14.7 


1026 


Allegan 


H5 


12.6 


1096 


Allegan 


H5 


12.2 


1003 


Nadiabondi 


H5 


11.7 


1087 


Nadiabondi 


H5 


12.0 


1001 


Forest Vale 


H4 


3.64 


897 


Forest Vale 


H4 


1.85 


813 


Ste. Marguerite 


H4 


0.638 


764 


Ste. Marguerite 


H4 


0.396 


713 



mechanism that may be responsible for the low value of sr f is that our optimized value represents a mean value of the actual 
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Fig. 5. Evolution history of the value of x 2 during the optimi- 
sation process with progressing number of generations, n, for 
the model where the heat conductivity at zero porosity, K\» is 
included in the optimisation process. Shown is for each genera- 
tion the particular individual with the lowest value of x 2 of all 
individuals of that generation, the formation time, ff orm , of this 
individual, and the radius of the body, R. For two representative 
meteorites, Kernouve and Allegan, the evolution of the optimum 
value of the burial depth, d, is also shown. 



value of <f> SI f over the initial heating and the extended subsequent 
cooling period of 100 Ma. During that period impacts certainly 
continue with a decreasing rate that do not contribute to a signif- 
icant further growth, but that form a regolith surface. According 
to m odel calculations this has a satu ration porosity of 20% to 
30% dRichardso rl2009tlWarrenll201 lb that would fit to the value 
of sr f = 0.20 found by the optimisation process. 

5.4. Model with optimised heat conductivity 

The pre-factor K\, in our approximation for the porosity- 
dependent heat conductivity, Eq. ([6j, corresponds to the heat 
conductivity of the bulk material (i.e., at <p — or D = 1), deter- 
mined by e xtrapolating data obtained f rom a couple of ordinary 
chondrites (Yomogida& Matsui 1983) to zero porosity. Figure 
6 in that paper and our Fig. [6] shows that there is considerable 
scatter in the data which results in some uncertainty of the value 
of K\,. The value for K\, used in the above model calculation 
(see Table [5]l may therefore be an unfavourable choice from the 



-e- 




Fig. 6. Data for thermal conductiv ity K exp of ordinary chondrites 
(from Yomogida & Matsui 1983) divided by the analytical fit 
K(<p) given by Eq. © using a value of K\, = 2.652 obtained as 
optimum value for fitting the cooling history. Full black circles 
represent data for H chondrites, open circles data for L chon- 
drites. The dashed line corresponds to the value of K = 4 used 
in model with fixed K. 



allowed range of values compatible with experimental values. 
Additionally, there is the problem that the small suite of chon- 
drites used for determining K\, may not be representative for the 
bulk of the material of the parent body. 

Therefore we performed a model calculation where K\, was 
considered to be a free parameter of the problem for which an 
optimal value is determined by the optimisation procedure. The 
remaining number of data-points for closure times with which to 
compare is still big enough to allow for an additional parameter, 
because we now have 28 data points and 14 parameters (6 for 
the body and the 8 burial depth's). 

The model, again, used 250 mass-shells. The optimisation 
was run for about 3 244 generations with hundred individuals 
per generation. The finally accepted optimum fit had a value of 
X 2 - 6.033 in generation 3229. With this value forx 2 the fit can 
also be considered as an excellent one in view of the fact that an 
acceptable fit would require only x 2 5= 14 (cf. Sect. 14.6b . How 
the value of x 2 evolved during the optimisation process for this 
model is shown in Fig. [5] The somewhat "noisy" appearance of 
the evolution history reflects the permanent mutations and cross- 
ing of genes in the evolution algorithm that not always results in 
children that are as fit as their parents. 

The model parameters of the final optimum fit for the prop- 
erties of the body and the burial depth's of the meteoritic mate- 
rial are given in Table [5] and the model results for the individual 
meteorites are shown in Fig. [3] The already good fit to the data 
points shown in Fig.|2]is once more slightly improved. The two 
meteorites Nadiabondi and Ste. Marguerite, however, drop out 
again and are not well fitted. 

The thickness of the remaining porous outer layer is also 
given in Table|6] This layer is defined as that region where <p(r) > 
0.1; its lower boarder is given by the radius where <p(r) = 0.1. 
Around this border the heat conductivity K(<p) drops nearly dis- 
continuously from the bulk value to the low value K((f> SI f). 
The blanketing effect of this layer is responsible for the shallow 
burial depths of the H4 and H5 chondrite s found in the model, 
which are also found in other models (e.g [Bennett & McSweenl 
I1996L lAkridge et alJll998t lHarrison & Grimmll2010l) that intro- 
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Fig. 7. Results for the individual meteorites if a Miyamoto-like model is used for the model fit. The burial depth's of the meteorites 
found for this model are given in Table[8] For meaning of symbols see caption to Fig. [2] 



Table 8. Parameters and burial depths of the meteorites for the 
Miyamoto-like optimised H chondrite parent body model. 



Quantity 


value 


unit 




Radius 


R 


143.67 


km 


Formation time 


f form 


2.157 


Ma 


Heat conductivity 


K b 


4.00 




60 Fe/ 56 Fe ratio 




2.07 x 10~ 7 




Surface temperature 




236.2 


K 


Porosity 


(/>srf 


8.0% 




Maximum temperature 


T c 


1 175 


K 


Porous outer layer 




none 


km 


Meteorite 


type 


depth [km] 


T'max [K] 


Estacado 


H6 


65.7 


1 150 


Guareiia 


H6 


63.5 


1 149 


Kernouve 


H6 


51.6 


1 136 


Richardton 


H5 


21.5 


1026 


Allegan 


H5 


20.1 


1013 


Nadiabondi 


H5 


18.8 


1003 


Forest Vale 


H4 


8.29 


828 


Ste. Marguerite 


H4 


3.98 


691 



duce an insulating outer regolith layer instead of modelling the 
sintering process. 

In this model the radius of the body is smaller than in the 
model with prescribed K\,. With R = 123km it is 30km bigger 



than the average radius of 93 km of (6) Hebe. If (6) Hebe really 
is the parent body, this would mean that material corresponding 
bf to 24% of the initial radius or 57% of the initial mass has 
been knocked-off by collisions. This model seems to be more 
plausible because the body then retained nearly half of its initial 
mass and did not lose most of its mass as in the preceding model. 
This makes it more plausible that we still find H4 chondrites. 

The value for K\, required for an optimum fit is significantly 
lower than the average value for K\, obtained from the labora- 
tory data but still inside the range of experimental values. This is 
shown in Fig. |6] wh ere we plot the ratio K exp /K((p) for all chon- 
drites from lYomogida & Matsuil d!983l) . If one would have used 
only the experimental values for the H chondrites, the optimised 
value for K\, would be definitely too low. This may result from 
the fact that, as already m entioned, the H chondrites used by 
Yomogida & Matsui| (Il983l) for measuring the heat conductivity 
and its dependence on <p are unfortunately not representative for 
the parent body as a whole. For instance could the material in 
the more central part of the body have for some reason differ- 
ent properties than the zones probed by the meteorites. Another 
probability could be that we neglected the temperature varia- 
tion of the heat conductivity. If one looks at Fig. [6] this variation 
seems to be small, but in fact the heat conductivity drops with in- 
creasing temperature and may be significantly lower at 1 000 K 
than the values seen in Fig. [6] Therefore the low value found in 
the optimisation may also simply mean that if we optimise K\,, 
the resulting value is kind of a "best average" between the lower 
conductivity of the hot core region of the parent body and of the 



12 



Stephan Henke et al.: Thermal history modeling of the H chondrite parent body 



higher conductivity in the cooler mantle region. The temperature 
dependence of the heat conductivity therefore should be consid- 
ered in future model calculations. 

5.5. Model with given Radius of (6) Hebe 

We also tried a model fit with the present-day radius R = 93 km 
of (6) Hebe to check whether a somewhat sub-optimal but, nev- 
ertheless, acceptable fit could be obtained for this case. The 
value of K\, is included in the optimisation process, but R is now 
fixed. The model, again, used 250 mass-shells. The optimisation 
was run for about 5 500 generations with hundred individuals 
per generation. The finally accepted optimum fit had a value of 
X 2 = 6.543 in generation 1 633. With this value for^ 2 the fit is 
not worse than the fit of the model with fixed R in view of the fact 
that an acceptable fit would require only^f 2 < 15 (cf. Sect. 14.61 ). 
The resulting set of parameters is given in Table [7] and the tem- 
perature variation at the burial depths of the meteorites is shown 
in Fig. |U 

The main difference between the models with optimised R 
and fixed R is that the resulting r sl f in the first model is 225.0 
and 306.5 K in the latter, all other things being very similar. 
The higher surface temperature is required to prevent a too rapid 
cooling of the outer layers of the smaller body. In practice this 
means that there is some ambiguity in the models with respect 
to the location of the body (surface temperature should decrease 
with increasing distance to the sun) and its radius. The optimi- 
sation process cannot really discriminate between the two cases 
considered here. 

5.6. Nadiabondi and Ste. Marguerite 

For Nadiabondi and Ste. Marguerite, some data points need to 
be discussed separately. 

For Ste. Marguerite the temperature increase during meta- 
morphism was insufficient for resetting some of the radioiso- 
topic systems, making it difficult t o constrain the cool ing his- 
tory of this meteorite. For instance, iKleine etaL I d20Q8h argued 
that the closure temperature of the Hf-W system is much higher 
than the metamorphic peak temperature of Ste. Marguerite. The 
Hf-W age of 1 .7 ± 0.7 Ma after CAI formation, therefore, must 
date a high-temperature event prior to the metamorphism. This 
event most likely was chondrule formation. This interpretation 
is consistent with our finding that the H chondrite parent body 
probably accreted at about 1.88 Ma after CAI formation. 

Second, some data of relatively poor precision do not really 
constrain the fit procedure, e.g., in the case of Nadiabondi the 
low temperature datum of fission track retention of the phos- 
phate merrillite. Nominally it yields an - impossible - negative 
cooling rate, and only within errors just relatively fast cooling 
can be inferred. Similarly, both fission track and Ar-Ar data of 
Ste. Marguerite have large uncertainties, precluding firm con- 
straints on the cooling history, although the nominal values plot 
on the fitted cooling-curves. 

5. 7. A Miyamoto-like model 

For comparison, the optimization was also done for a model 
where the porosity was set to a value of <p ST f = 0.08. This 
is a typical value as it is observed for various H6 chon- 
drites that have not been subject to strong shock mo difications 
dYomogida & Matsuill98llConsolmagno "et al.l2008l) . With this 
initial porosity the body is practically already compacted from 



the beginning on and sintering does not play any role. The heat 
conductivity then is practically constant throughout the body and 
for all times. The prefactor K\, of the heat conductivity was set 
to the same value as in Sect. l5.3l and was not included in the op- 
timisation process. Such a model corre sponds essentially to th e 
kind of model onsidered originally by iMivamoto et alj d 198 lb . 
the mayor difference being that heating by 60 Fe is included. 

The model, again, used 250 mass-shells. The optimisation 
was run for about 3 260 generations with hundred individuals per 
generation. In generation 2299 generations an optimum value 
of x 2 = 5.94 was obtained that was considered as sufficiently 
good. The model parameters of the optimised model are shown 
in Tableland Fig.Qshows the fit for the individual meteorites. 

It is obviously possible to obtain also an acceptable fit to 
the observed cooling history of the meteorites with such kind of 
model. The parameters of such a model are somewhat different 
from the preceding cases because of the lack of an insulating 
porous outer layer and a higher heat conductivity during the ini- 
tial warm-up phase. A significant contribution of 60 Fe decay to 
the heating is required to compensate for the more easy heat-loss 
in the present case. Though the model as a whole does not look 
worser than the model with optimised K^, it is physically not 
acceptable since the lack of a porous outer layer is unrealistic. 

5.8. Comparison with earlier results 

There are a number of previous attempts to reconstruct the prop- 
erties of the H chondrite parent body using the empirically de- 
termined cooling rates of H chondrites derived from different 
depth within the parent body. Table [9] gives an overview on the 
results and the values of some key parameters of some of these 
models. The table lists all the models that attempted to obtain in 
some way a "best fit" between properties and cooling histories 
derived from laboratory studies and a thermal evolution model 
of the parent body. The results are generally the putative radius 
and formation time of the body, and the burial depths of the me- 
teorites. All studies, including ours, showed, that the observed 
properties of the H chondrites are compatible with a parent body 
of about 100 km radius and an "onion shell" model. In this model 
the body is first heated within a period of ca. 3 Ma by decay of 
radioactive nuclei, mainly 2f Al, and then cools down over an ex- 
tended period of ca. 100 Ma, without suffering catastrophic col- 
lisions that would have disrupted the body, thereby disturbing its 
thermal evolution record. 

Despite the different assumptions made in the model calcu- 
lations, from the rather simplistic model of a homogeneous body 
with constant material p roperties, that can be treated analytically 
(Miyamot o et al.lQ98lh . to models with rather comp lex imple- 
mented physics like that of lHarrison & Grimml J20 1 0) and ours, 
the derived properties of the parent body are similar. These do 
not seem to depend very critical on the particular assumptions. 
The more strongly deviating results for the burial depths of the 
meteorites cannot be used to discriminate between models, be- 
cause no independent information on this is presently available. 
This would require that methods are developed to determine the 
maximum pressure experienced by a meteorite. 

The most important difference between our model and previ- 
ous models is that we included sintering of the granular material 
in the modelli ng. This process is not cons idere d in other mod- 
els, ex cept for lYomogida & Matsuil (1 1984h and|Benoit & Searsl 
(119961) . The model oflYomogida & Matsui ( 1984) unfortunately 
does not consider the important contribution of 26 A1 to the heat- 
ing. lTierefore_tlris model cannot be compared to other mod- 
els. iBenoit & Searsl (Il996h gave only very few details on their 
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Table 9. Properties of the parent body of H chondrites derived by different thermal evolution models, and some key data used for the 
modelling. In some models an insulating outer regolith layer is assumed; in that cases a core radius and layer thickness are given. 



Model 


R 

km 


Ma 


K 


M 
kg 


K 
W/mK 


Cp 

J/kgK 


p 

kg/m 3 


h f 
W/kg 




Method 


Sintering 


Miyamoto et al. (19811 


85 


Z.J 




O.Z 1U 


i 
i 


OZj 


J zuu 


1 .OZ 1U 


-7 


analytic 




Benoit& Sears (1996) 


99.0+1.0 


2.2 


200 


n.a. e 


variable d 


variable b 


3 780 a 


n.a. e 




numeric 


yj + regolith 


Akridge et al. Q998) 


97.5+2.5 


2.3 


200 


1.4 10 19 


variable 11 


variable c 


3 450 






numeric 


regolith 


Trieloff etal. (2003) 


100 


2.5 


300 


1.3 10" 


1 


625 


3 200 


1.82 10~ 


-7 


analytic 




Kleine et al. (2008) 


100 


2.7 


300 


8.2 10 18 


1 


625 


3 200 


1.8210" 


-7 


analytic 




Harrison & Grimm (2010) 


99.2+0.8 


2.2 


170 


1.4 10 19 


variable 11 


variable b 


3 250 


1.8210" 


-7 


numeric 


regolith 


This work 


123 


1.88 


225 


3.0 10 19 


2.65 a 


variable b 


3 780 a 


1.67 10 


-7 


numeric 


V 


(6) Hebe g (today) 93 

(diameter 205 x 185 x 170 km) 






1.28 10 19 
±6.4 10 17 






3810 
±260 











Notes: (a) for consolidated material, varies with porosity (b) lit of c p to meteoritic data from Yomogida & Matsui 1 1984), at T = 500 the value is 865 J/kgK (c) forsterite, (d) fit of K to 
meteoritic data from Yomogida & Matsui 1 1984), at T - 500 the value of the consolidated material is 3.66 W/mK and varies only slightly with temperature, (e) no or insufficient data 
presented, (f) at time of CAI formation, given value for 26 A1 only, some models consider also the less efficient heat sources ( fi0 Fe, long lived radioactives), (g) from data collection Baer et al. 
<20Tll) . 



model results such that one cannot make a compa rison with our 
and ot her mod els also in this case. Th e models of Akridg e et al.l 
(1998)and|H arrison & Grimrnl dxHHOl) do not consider sintering, 
but simulate the insulating effect of a porous outer layer, that 
escaped sintering, over a consolidated body by introducing two 
outer layers with low filling factor, a "regolith" and a "mega- 
regolifh" layer. The thickness's of both layers were assumed by 
Akri dge et all (|l_998b but they were dete r mined by an optimisa- 
tion of the model bv lHarrison & Grimrnl J2010) and a thickness 
of 410 m was obtained for each of the layers. This compares well 
with a thickness of < 1 km of a residual porous outer layer in our 
model. The thermal shielding of the consolidated interior of the 
body by a thin insulating outer layer of low thermal conductivity 
is obviously essential for obtaining a good fit between thermal 
ev olution models to meteoriti c data as it is obtained in the model 
bv lHarrison & Grimrnl d20 10) and in ours. 

A major difference between our model results and others is 
the earlier derived formation time, ff olm . This difference probably 
results from the lower mass-density used in other model calcu- 
lations and the resulting lower mass of the body. Because of the 
resulting higher total heat capacity of the whole body ( Jpc v dV) 
in our model more heat from decay of radioactive nuclei is re- 
quired to heat the body to the temperatures as witnessed by the 
observed cooling history of the various meteorites. Hence, an 
earlier formation time is required for our optimised model since 
the specific heating rate at zero time (equal to time of CAI for- 
mation) is almost the same in all models. The higher average 
density in our model is a result of the sintering process which 
results in very low or vanishing porosity over most of the vol- 
ume of the body. 

If (6) Hebe really is the parent body, the mass-density in our 
model fits well to that determined for Hebe (I James & Cheslevl 
20081 lBaeretld1l2011h . The densities used in other model cal- 
culations then would be too low and the earlier formation time 
found in the present model calculation therefore is probably 
more realistic than the later formation times found in other mod- 
els. 

The surface temperature, r sr f , is determined in our model by 
the optimisation process, in the other model it seems to be cho- 
sen more or less arbitrary within the limit expected for bodies in 
the asteroid belt and in some cases adjusted to improve the model 
fit between closing temperatures and closing times observed for 
meteorites. The dependence of the thermal evolution model on 



the value of r sr f is generally weak (e.g. lHenke et al.ll2012l) . The 
value determined by our optimisation procedure hinges mainly 
upon the fit to the two H4 chondrites that are most strongly af- 
fected by r sr f because of the shallow burial depths required by 
their petrologic type. 

There are a n umber of other model calculations 
(Bennett & McSweenl 



Sahii oal et alJ 120071: |W: 



19961 



Hevev & Sanders! l2006t 
1) that also studied the 



thermal evolution of asteroids and give important insight into 
their properties and the processes operating in them, but did 
not attempt to obtain definite results for the properties of the 
parent body of H chondrites from some kind of "best fit" to 
meteoritic data. Only the range of possible values is studied in 
these papers. This does not enable a direct comparison, but the 
range of possible values for radii and formation times given are 
in accord with the results shown in Table [9] 

5.9. Orbit of parent body 

If the number obtained for r sr f from our optimisation process is 
taken seriously, it measures the average surface temperature for 
the first = 20 ... 30 Ma of evolution where the temperature for 
Forest Vale and Ste. Marguerite was above the closure tempera- 
ture for Pu fission tracks. Since the observed lifetime of accretion 
disks is less than 10 Ma (Hai sch et al.l l2()01) the surface temper- 
ature is determined by illumination by the young proto-sun and 
thermal emission from the surface, in particular during the most 
critical phase where the temperature at the burial depth evolved 
through the closing temperature for Pu fission tracks. Then we 
can estimate the average radius of the orbit of the parent body 
from the relation (for a rapidly rotating body) 



(l-A)S(l) 
4a 2 



(13) 



where cr is the Stefan-Boltzmann constant, A the albedo of the 
surface, a the distance to the protosun in AU, and S(l) the solar 
constant at 1 AU at this early evolution period. 

For the solar constant S(l) we have to observe that the sun 
ignites at about 10 Ma hydrogen burning and has the luminos- 
ity of a 1 M0 star on the zero age main sequence. Its luminos- 
ity during t he first 100 Ma is by about 30% less than its present 
value (e.g. iD'Antonia & Mazziteililll994l) . The typical surface 
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albedo for S asteroids is 0.2 ( Rvan ~&" Woodwardl20 1 Oh . The sur- 
face albedo of (6) Hebe i s particular high with 0.291 ± 0.029 
dRvan & Woodwardll2010h . The albedo found today is probably 
not typical for the material that covered the surfaces of the bodies 
4.5 Ga ago, because highly processed material and the products 
of its grinding are presented at the surface. It must had contained 
a lot of more fine-grained material that only later matured by 
crystal growth at elevated temperatures and probably was darker. 

We assume arbitrarily A = 0.1 and find a radius of the orbit 
of a = 1.23 AU (an albedo of A = 0.2 yields a = 1.16 AU). This 
is smaller than the present average radius of the orbit of 2.4 AU. 
If the temperature r sr f found in optimisation is really represen- 
tative for the average surface temperature during the essential 
phase of the thermal evolution, and if (6) Hebe is the parent body 
of the H chondrites, then our result means that (6) Hebe formed 
closer to sun than its present distance. However, a more detailed 
treatment of the processes determining the surface temperature 
is necessary, before any definite conclusions can be drawn. 

6. Concluding remarks 

We combined the code for modelling the internal structure and 
evolu tion of planetesimals of the 100 km size class dHenke et al.l 
l20T2h with a code for optimisation of a quality function that is 
based on an evolution algorithm. This allows finding a consistent 
set of parameters that describe the properties and evolution of a 
planetesimal and the burial depths of the meteorites. A such the 
thermal evolution of the model body reproduces very closely the 
empirically determined cooling history of all H chondrites with 
sufficient available data. 

The results show, that the "onion shell" model successfully 
reproduces the thermal evolution history of all the meteorites 
used in this study. For six of the eight meteorites a burial depth 
can be found for which the empirical cooling data can be fitted in 
an excellent way with the local evolution of temperature at that 
burial depth. For the remaining two cases it is still possible to 
find a burial depth that results in a good fit if the Pu fission-track 
datum of Nadiabondi is rejected. This confirms earlier conclu- 
sions that the thermal evolution of the H chondrite parent body 
and its initial thermal structure has not been disturbed by catas- 
trophic impacts. 

The reconstructed properties of the parent body of the H 
chondrites are comparable to what has been found in other model 
calculations. The different published models treated the problem 
in different kinds of approximations, solving the heat conduc- 
tion equation being almost the only feature common to all mod- 
els, and used different values for the material properties. Our 
estimated size of 123 km of the parent body is somewhat bigger 
compared to the 100 km found in other models, and requires a 
slightly earlier formation time of 1.88 Ma compared to the 2.2 - 
2.5 Ma found in other models. Such minor differences are cer- 
tainly due to the different approaches used. 

Our study shows, that the input physics of the models has to 
be improved, before more reliable data on the parent body can 
be obtained. Such improvements are currently underway in our 
group. In particular the heat conductivity has to be calculated 
from the specific properties of chondritic material. It will also 
be necessary to obtain more precise radiometric ages for a wider 
range of samples. 

The advantage of our method is that it allows for a simul- 
taneous and consistent fit of several model parameters to a big 
number of observational restrictions and is very flexible if the 
kind of data to be included in the optimization process is to be 
changed. A minor shortcoming of the optimisation method used 



is that a considerable number of complete evolution models has 
to be calculated during the course of the optimisation process, 
but our implicit solution method is fast enough that no restric- 
tion as to the applicability of the method to our problem results 
from this. 

Acknowledgements. We thank the referee S. Sahijpal for a constructive and 
helpful referee report. This work was supported in part by 'Forschergruppe 
759' and in pail by 'Schwerpunktprogramm 1385'. Both are supported by the 
'Deutsche Forschungsgemeinschaft (DFG)'. S.H. is Member of the IMPRS for 
Astronomy & Cosmic Physics at the University of Heidelberg. 



References 

Akridge. G.. Benoit. P. H., & Sears, D. W. 1998, Icarus, 132, 185 

Alexander, C. M. O. 2005, in Astron. Soc. of the Pacific Conf. Ser., Vol. 341, 

Chondrites and the Protoplanetary Disk, ed. A. N. Krot, E. R. D. Scott, & 

B. Reipurth, 972-1002 
Arzt, E., Ashby, M. F, & Easterling, K. E. 1983, Metallurgical Transact. A, 14A, 

21 1 

Baer, J., Chesley, S., & Britt, D. 201 1, Asteroid Masses V2.0., eAR-A-COMPIL- 

5-ASTMASS-V2.0, NASA Planetary Data System 
Bennett, M. & McSween, H. Y. 1996, Meteoritics & Plan. Sci., 31, 783 
Benoit, P. & Sears, W. G. 1996, Lunar Planet. Sci. Conf. Lett., 28, 1178 
Birck, J. L. & Lugmair, G. W. 1988, Earth & Plan. Sci. Lett., 90, 131 
Bouvier, A., Blichert-Toft, J., Moynier, F, Vervoort, J., & Albarede, F. 2007, 

Geochim. Cosmochim. Acta, 71, 1583 
Charbonneau, P. 1995, ApJS, 101 

Cherniak, D., Lanford, W. A., & Ryerson, F. J. 1991, 

Geochim. Cosmochim. Acta, 55, 1663 
Consolmagno, G., Britt, D., & Macke, R. 2008, Chemie der Erde / Geochemistry, 

68, 1 

D'Antonia & Mazzitelli. 1994, ApJS, 90, 467 

Dauphas, N., Cook, D. L.. Sacarabany, A., et al. 2008, ApJ, 686, 560 
Dodd, R. 1969, Geochim. Cosmochim. Acta, 33, 161 
Dodson, M. Ft. 1973, Contributions to Mineralogy and Petrology, 40. 259 
Fei, Y, Berfka, C. M., & Finger, L. W. 1997, Science, 275, 1621 
Fischmeister, H. F. & Arzt, E. 1983, Powder Metallurgy, 26, 82 
Gaffey, M. J. & Gilbert, S. L. 1998, Meteoritics and Planetary Science, 33, 1281 
Ganguly, J. & Tirone, M. 2001, Meteoritics and Planetary Science, 36, 167 
Gopel, C, Manhes, G., & Allegre, C. J. 1994, Earth & Plan. Sci. Lett., 121, 153 
Guttler, C, Krause, M., Geretshauser, R. J., Speith, R., & Blum, J. 2009, ApJ, 
701, 130 

Haisch, Jr., K. E., Lada, E. A., & Lada, C. J. 2001, ApJ, 553, L153 
Harrison, K. P. & Grimm, R. E. 2010, Geochim. Cosmochim. Acta, 74, 5410 
Henke, S.. Gail, H.-P, Trieloff, M., Schwarz, W. H., & Kleine, T. 2012, A&A, 
537, A45, (paper I) 

Herpfer, M., Larimer, J., & Goldstein, J. 1994, Geochim. Cosmochim. Acta, 58, 
1353 

Hevey. P. J. & Sanders, I. S. 2006, Meteoritics & Plan. Sci., 41, 95 

Hopfe, W. & Goldstein, J. 2001, Meteoritics & Plan. Sci., 36, 135 

James, B. & Chesley, S. R. 2008, Celestial Mech Dyn Astr, 100, 27 

Kakar, A. K. & Chaklader, A. C. D. 1967, J. Appl. Phys., 38, 3223 

Kita, N., Nagahara, H., Togashi, S., & Morishita, Y. 2000, 

Geochim. Cosmochim. Acta, 64, 3913 
Kleine, T., Hans, U., Irving, A. J., & B., B. 2012, Geochim. Cosmochim. Acta 
Kleine, T., Touboul, M., Bourdon, B., et al. 2009, Geochim. Cosmochim. Acta, 

73,5150 

Kleine, T., Touboul, M., Van Orman, J., et al. 2008, Earth & Plan. Sci. Lett., 270, 
106 

Krause, M., Blum, J., Skorov. Y, & Trieloff, M. 201 la, Icarus, 214, 286 
Krause, M., Henke, S., Gail, H.-P, et al. 2011b, Lunar Planet. Sci. Conf. Lett., 
42, 2696 

Larsson, P. L., Biwa, S., & Storakers, B. 1996, Acta mater., 44, 3655 
McSween. H. & Patchen, A. 1989, Meteoritics, 24, 219 

Minster, J. F. & Allegre, C. J. 1979, Earth and Planetary Science Letters, 42, 333 
Miyamoto, M., Fujii, N., & Takeda, H. 1981, Lunar Planet. Sci. Conf. Lett., 12B, 
1145 

Morbidelli, A., Gonczi, R., Froeschlee, C, & Farinella, P. 1994, A&A, 282, 955 
Nagasawa, M., Thommes, E. W., Kenyon, S. J., Bromley, B. C, & Lin, D. N. C. 

2007, in Protostars and Planets V, ed. B. Reipurt, D. Jewitt, & K. Keil 

(Tucson: University of Arizona Press), 639-654 
Nyquist, L. E., Kleine, T., Shih, C.-Y, & Reese, Y. D. 2009, Geochim. 

Cosmochim. Acta, 73, 5115 
Pellas, P., Fieni, C, Trieloff, M., & Jessberger, E. 1997, Geochim. Cosmochim. 

Acta, 61, 3477 



15 



Stephan Henke et al.: Thermal history modeling of the H chondrite parent body 

Quitte, G., Markowski, A., Latkoczy, C, Gabriel, A., & Pack, A. 2010, ApJ, 720, 
1215 

Rao, A. S. & Chaklader, A. C. D. 1972, J. Americ. Ceramic Soc, 55, 596 
Rechenberg, I. 1973, Evolutionsstrategie - Optimierung technischer Systeme 

nach Prinzipien biologischer Evolution (Stuttgart: frommann-holzboog) 
Rechenberg, I. 1989, in Machine Learning: Principles and Techniques, ed. 

R. Forsyth, (London: Chapman and Hall), 83-103 
Richardson, J. E. 2009, Icarus, 204, 697 

Rudraswami, N. G. & Goswami, J. N. 2007, Earth & Plan. Sci. Lett., 257, 231 
Rudraswami, N. G, Goswami, J. N., Chattopadhyay, B., Sengupta, S. K., & 

Thapliyal, A. P. 2008, Earth and Planetary Science Letters, 274, 93 
Rugel, G, Faestermann, T., Knie, K., et al. 2009, Phys. Rev. Lett., 103, 072502 
Ryan, E. L. & Woodward, C. E. 2010, AJ, 140, 933 

Sahijpal, S., Soni, P., & Gupta, G. 2007, Meteoritics & Plan. Sci., 42, 1529 
Slater-Reynolds, V. & McSween, J. 2005, Meteoritics & Plan. Sci., 40, 745 
Storakers, B., Fleck, N. A., & McMeeking, R. M. 1999, J. Mech. Phys. Solids, 
47, 785 

Tang, H. & Dauphas, N. 2012, Lunar Planet. Sci. Conf. Lett., 43, 1703 
Taylor, G. J., Maggiore, P., Scott, E. R. D., & Keil, A. E. 1987, Icarus, 69, 1 
Torppa, J., Kaasalainen, M., Michalowski, T, et al. 2003, Icarus, 164, 346 
Trieloff, M., Jessberger, E. K., Herrwerth, I., et al. 2003, Nature, 422, 502 
Warren, P. H. 2011, Meteoritics & Plan. Sci., 46, 53 

Weidenschilling, S. J. & Cuzzi, J. N. 2006, in Meteorites and the Early Solar 

System II, ed. D. S. Lauretta & H. Y. McSween (Tucson: University of 

Arizona Press), 473^185 
Wittmann, A., Swindle, T. D., Cheek, L. C, Frank, E. A., & Kring, D. A. 2010, 

J. Geophys. Res., 1 15, E07009 
Yomogida, K. & Matsui, T. 1983, J. Geophys. Res., 88, 9513 
Yomogida, K. & Matsui, T. 1984, Earth & Plan. Sci. L., 68, 34 



16 



